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In these lectures, I describe the techniques used within the QCD sum rule ap- 
proach. The basic concepts of the approach are introduced using a simple model 
of quantum-mechanical oscillator in 2+1 dimensions. Then I discuss their field- 
theoretical extension and the construction of the operator product expansion for 
current correlators in QCD. The calculation of static parameters is illustrated on 
the example of sum rules in the vector and axial meson channels. Finally, the 
QCD sum rule calculation of the pion electromagnetic form factor is presented as 
an example of application of the method to dynamic characteristics of the hadrons. 



1 Introduction 

1.1 Introductory remarks 

QCD sum rnlfls, invented more than 20 years ago by Shifman, Vainshtein 
and Zakharo\o are still a rather popular approach for theoretical study of the 
hadronic structure. Applied originally to the simplest static hadronic char- 
acteristics, like masses, leptonic widths etc., they were also used to calculate 
much more complicated things like hadronic wave functions and form factors. 
To begin with, I will briefly discuss the basic questions 

• Why we need QCD sum rules? 

• What are QCD sum rules? 

• How the QCD sum rules are working? 

using for illustration the QCD calculation of the pion form factor. In subse- 
quent sections, I will discuss in more detail the machinery of QCD sum rules in 
several cases: calculation of the properties of the lowest resonances in p-meson 
and 7r/Ai channels, and analysis of the pion form factor in various kinematical 
situations. 



1.2 Pion form factor in QCD 

According to the factorization theoreni'i, for sufficiently large momentum 
transfers, one can represent the pion form factor as a sum of terms of increasing 
complexity (see Fig.l). The first term (purely soft contribution, Fig.^a con- 
tains no short-distance (SD) subprocesses. The second term Fig.|l|b contains a 
hard gluon exchange. There are also corrections to the hard term: higher or- 
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Figure 1: Structure of factorization for the pion form factor in QCD 

der corrections Figjl]; containing extra as factors and higher twist corrections 
Fig.|l|d. In perturbative QCD one can calculate only the hard terms of this 
expansion. 

The first result of perturbative QCD for exclusive processes is the pre- 
' n for the asymptotic behavior of the pion electromagnetic form factor 

F^{Q^) ^ [ dx f t{xP, (1 - x)P; yP', (1 - y)P' , q) ^{x) ^{y) dy , (1) 
Jo Jo 

where f{x) is the pion wave function giving the probability amplitude to find 
the pion in a state where quarks carry fractions xP and (1 — x)P of its lon- 
gitudinal momentum P; q = P' — P is the momentum transfer to the pion 
and T(. . .) is the perturbatively calculable short-distance amplitude. The wave 
function 1^9(2;) accumulates long-distance information about the pion structure 
and cannot be calculated within the pQCD approach. Normally, some phe- 
nomenological assumptions are used about the form of this function depending 
on a single parameter x. Still, there is a challenge to calculate it from first 
principles of QCD. 

To calculate the soft contribution to the form factor one should know 
something like the light-cone wave function ^p{x,kj_) depending also on the 
transverse momentum k±_ . I say "something like" , because in principle it is 
impossible to prove that the soft term reduces in QCD to the convolution 

F(g2) ^ J ^p{x, kx)^*p{x, k± + xq) d^k±dx. (2) 

In fact, one should treat the soft term as a nonfactorizable amplitude and the 
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problem is to calculate it in some reliable way. 

The QCD sum rules are just giving the method to calculate the nonpertur- 
bative hadronic characteristics incorporating the asymptotic freedom property 
of QCD in a way similar to pQCD. 

1.3 Basic ideas of the QCD sum rule approach 

Among the existing approaches to the analysis of the nonperturbative effects 
in QCD the most close to perturbative QCD is the QCD sum rule method 
EJ. Let us formulate its basic ideas within the context of the pion form factor 
problem. 

It is evident that one cannot directly study the soft contribution with the 
on-shell pions, because then only long distances are involved. But perturbative 
QCD can be applied in a situation when all relevant momenta pi, p2 are 
spacelike and sufficiently large: \q^\, \p\\, > 1 GeV^. To describe the virtual 
pions one should use some interpolating field, the usual (for the QCD sum rule 
practitioners) choice being the axial current = d757"w. Its projection onto 
the pion state |P, tt) is proportional to /^r : 

(0|j5"|P,7r) ^z/^P". (3) 

Via the dispersion relation 

T{pi,pi,q^) = — dsi ds2 -, 4 

Jo Jo (si - Pl)(S2 ~ P2) 

one can relate the amplitude T{p1,p2, q^) to its time-like counterpart p{si, S2j 9^) 
containing the double pole term 

P.(si,S2,g') = 7T^f^S{si ~ ml)S{s2-ml)F^{Q^) (5) 

corresponding to the pion form factor. However, the axial current has nonzero 
projections onto other hadronic states (Ai— meson, say) as well, and the spec- 
tral density p{si,S2,q^) contains also the part /9''*^'''^'''**''*^*(si, S2, g^) related 
to other elastic and transition form factors. This is the price for going off the 
pion mass shell. The problem now is to pick out the F^^ term from the whole 
mess. 

Of course, calculating T{p\,p\,q^) in the lowest orders of perturbation 
theory one never observes something like the pion pole: one obtains a smooth 
function p^'^^^isi, S2,q^) corresponding to transitions between the free-quark 
Mc?— states with invariant masses si and S2, respectively. The difference be- 
tween "exact" density p(si, S2,q'^) and its perturbative analog pP'^^*{si, S2,q^) 
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Figure 2: Structure of the operator product expansion for T(pl,pl,q'^) 

is reflected by additional nonperturbative contributions to T{p\,p^, q^). These 
contributions are due to quark and gluon condensates (gg), (GG) etc., de- 
scribing (and/or parameterizing) the nontrivial structure of the QCD vacuum 
state. Formally, these terms appear from the operator product expansion for 
the amplitude T{pl,pl,q^) (Figj|): 

T{plplq^)=Tr'^^\plplq^)+a^^ + b^^ + ... . (6) 

The problem now is to construct such a model of the spectral density 
p{sii S2,q^) which gives the best agreement between two expressions for T, i.e., 
Eqs . and . Naturally, having only a few first terms of the 1 /p^ — expansion 
one can hope to reproduce only the gross features of the hadronic spectrum in 
the relevant channel. Still, just using the simple fact that the condensate con- 
tributions die out for large p^ , one obtains the global duality relation between 
quark and hadronic densities 

/•OO 

dsi dS2{p{si,S2,q^)~pP''^'{si,S2,q^))=0. (7) 

Jo 

In other words, integrally the two densities are rather similar. One can even 
expect that they are very close for high s to secure the convergence of the 
integral (0). 

The standard ansatz is to approximate the higher states contribution into 
p{si,S2,q^) by the free quark density (compare to the quasiclassical approxi- 
mation hor high levels in quantum mechanics): 

pMgher states ^^^^^^^^2^ = [1 - < So)0(s2 < Sq)] Z^''* (si , S2 , ) , (8) 

with So being the effective threshold for the higher states production. 



4 



1.4 Summary 

To summarize, the basic ingredients of the QCD sum rule approach are 

• Operator product expansion, or factorization of the short-distance, perturba- 
tively calculable amplitudes, with the long-distance information accumulated 
by the quark and gluon condensates - the basic parameters of the approach - 
characterizing the properties of the QCD vacuum. 

• Dispersion relations, relating the perturbatively calculated amplitude with 
the physical hadronic parameters. 

• To extract the hadronic parameters one should construct an ansatz for the 
hadronic spectrum (like "first resonance plus continuum", etc.) and then 
perform fitting procedure. 

Just from the above list it is clear that QCD sum rules cannot produce 
exact results. Their precision is unavoidably limited by the following facts : 

• One can calculate normally only few first terms of the OPE. 

• The model ansatz reproduces the real density only approximately. 

• Furthermore, in some channels could appear nonperturbative short-distance 
contributions ("direct instantons" ) . 

As a result, using the QCD sum rules, one can normally expect only 10% 
- 20% accuracy. 

2 QCD Sum Rules in Quantum Mechanics 

We are going to apply the QCD sum rule method in the theory, where no 
exact results for hadrons were obtained so far. But, as formulated above, the 
QCD sum rule method is quite general. In fact, many features of the method 
can be demonstrated on the exactly solvable model of the quantum-mechaniC|al 
oscillator. In the discussion below, -we basically follow the pioneering papefl, 
with some simplifying modificationicro. 

2.1 Sum rules for the lowest state of two-dimensional oscillator. 

Let us consider the simplest example of a system with confinement, the oscil- 
lator, with the confining potential being 

Vi^l = ^r^- (9) 

Formulas are the simplest (without square roots, complicated factorial terms, 
etc.) if we choose the oscillator in two spatial dimensions (-1- time dimension). 
In this case the energy levels are given by 

En = (2n + l)uj, (10) 
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and the values of the wave functions at the origin (quantum-mechanical ana- 
logues of Jtt and other leptonic decay parameters) are the same for all levels 
corresponding to radial excitations with zero angular momentum: 

\M0)\' = — . (11) 

TT 

To calculate Eq and "00(0) using the sum rule approach one should proceed 
in a rather peculiar way: instead of concentrating on the lowest level, one 
should consider the weighted sum of all levels: 

oo 

M{e)^J2 IV'.(0)pe-^^/^ (12) 

fc=0 

in which the lowest level dominates only for very small values of e. For large e 
all levels are essential. This combination (|l^) is directly related to the two-time 
Green function 

oo 

G(xi,ti|x2,t2) =^V'^(^2)Vfc(^i)e-'''''<*^"*^) (13) 

fe=0 

describing the probability amplitude for transition from the point Xi at time 
ti to the point X2 at time ^2- To get M{e) one should take a;2 = xi = 0, ii = 
and ^2 = l/ie. In other words, the transition is taken in the imaginary time. 

In our case, explicit form of Ajf(e) can be easily calculated using Eqs. (10) 
and (0): 

M-(.) ^ y e-(2^+i)-A ^ — (14) 

^— ^ V TT / 7rsinh(w/e) 

In the large- e limit one has 



M°'"'(e) = — 1 + + ...] (15) 

^ ' 2tt \ 6e2 360 15120 J ^ ' 

i.e., the first term of the large-e (=small-time) expansion has no dependence 
on uj, the paramater characterizing the strength of oscillator potential. This 
simply means that at small time differences the Green function coincides with 
the free one: 

and 

Mfree^^^ = ^ (17) 
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The corrections to the free-particle behavior in Eq. (^5|) are in powers of llP' . 
This corresponds to the perturbation expansion in the potential V{r) (^): 

M = M^'''='= + M^'''=^0F®M^'''='= + M^'''=^0F®M^'''=^0F(8)M^'''='= + . . . (18) 

The function Af/''^'=(e) can also be represented as a sum over the (free) 
states with — k"^ /2ra: 

M'-{e) = / ^e-^^/^ = f He-^'^dE^^- H e-^/^piE)dE . 

(19) 

We introduced here a new function 

pf^--(E)^'^e{E>Q), (20) 

the spectral density corresponding to the free motion, which will play an im- 
portant role in the subsequent discussion. The oscillator function M°'"^(e) also 
can be written in the spectral representation with the spectral density p°^'^{E) 
which is in this case a superposition of the delta-functions: 

oo 

p°'\E)=muj^5{E^{2k + l)uj). (21) 

At first sight, the two densities are completely different. However, if one in- 
tegrates them from zero to the midpoint between the first and second levels, 
one obtains the same result: 



-1 t^2u) -i f'2u) 

— / p^''^''(E)dE = uj=— p°"'{E)dE . 
m Jo m Jo 



(22) 



Moreover, similar relations hold for other levels: 

/ pf'-'^E)dE = muj = p°'^{E)dE . (23) 

J2kuj J2kL0 



The equality is not violated even if one multiplies p{E) by an extra power of 
the energy E: 

p2{k+l)w ^2(fe+l)w 

/ pf''''^{E)EdE ^muj^{2k+\) = I p°'"'{E)EdE. (24) 

J2kuj J2kuj 

Thus, there exists the local duality between each resonance and the free states. 
It looks like each level "takes" its part of the free spectral density, leaving the 
integral unchanged. 



Turning back to the cj-expansion for A/°'"^(e) (Eq. (15)) 

^ ' ^ ' \ 6e^ 360 £4 15120 e< 

one can write 

M°'"=(e) - M^^'^'^ie) ^ 0{Lu^/e) + Oioj^/e^) + . . . . (26) 
or, in terms of the spectral densities 

1 /-OO A 



f A 

/ [p-(iJ)-p/-(iJ)]e-^/^d£; = 5^4. (27) 

•^0 ^ 



The terms on the rhs, suppressed for large e by powers of 1/e are normally 
referred to as "power corrections" . Thus, the difference between the integrated 
spectral densities vanishes as e — > cx): 



1 f°° 
Jo 



(28) 



since there is no 0(1) term in the rhs of Eq. (p7|). Eq. ( p8| ) is a standard global 
duality relation. In quantum mechanics one can obtain such a relation for any 
potential regular at the origin. 

The lesson is that the two densities, despite the apparent absence of sim- 
ilarity, are rather close to each other in an integral sense. Now, the idea is to 
incorporate this similarity in the sum rule 

|^o(0)|^e-oA + ",.,,e. state,' ^ ^ ' ^ + Wol^ ' J^UOl^ + ' 

(29) 

by substituting the "higher states" contribution by that of the free states hav- 
ing energy above some threshold sq. The sum rule then reads 

\MO)\'e-^o/^ ^ i (i^^-so/A_^ + J_^_J]_^ + ^^^ (30) 
'^"^ 2V J 12e 720 £3 30240 ^ ' 

where 

|^o(0)P = -|V^o(0)|^ 
m 

The parameters to be extracted from the sum rule are Eq, V'ojO) and sq. The 
first step is to extract Eq. To this end we differentiate Eq. (^) with respect 
to (—1/e) to get 



2 .2 



\MO)\'Eoe-^"/^^'--'-{l+'-f) 



-sn/e ^"^ 7 LO^ 31 

e ^01'^ A 1 



12 240 e2 6048 

(31) 
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The ratio of the two equations just gives the expression for Eq: 



p ^ 2 2\'^^t'^ ' 12 240 7^ ' 6048 7^ ^ ■ ■ ■ (or)\ 

2V-^ ^ 7 12£ ^ 720 7^ 30240 ^ 

If we would take the whole series in uj/t and the exact expression for the 
higher states contribution, then the result for Eq would be the e-independent 
exact value i?o — lo. However, under approximations we made, the i?o-value 



extracted from Eq. (32) depends both on the auxiliary parameter e and the 
threshold parameter sq- The standard strategy is to take the value of sq, 
for which Eo has the minimal dependence on e in the intermediate region 
\ < e/uj < 2 where both the higher power corrections and the error due to 
the rough model of the higher states contribution are small (« 5 — 10%). 
Truncating the series by terms of the second order in uj^, one obtain^ from 
such a fitting sq = 1.6u; and Eq = 0.9w. Adding the third order terms improves 
the result: sq = l-75a; and Eq — 0.95w. In our case one can substitute into the 
sum rule the all-order result. The remaining e-dependence is then entirely due 
to the use of an approximate model for the higher states. For sufficiently small 
values of e all the curves give the exact value Eq = 2lo. The flattest curves are 
those with sq ^ 2uj. 

The second step is to extract the V'(O) value from the original sum rule 

(§) 



« e^°/^ ^ -(1 - e-^"/^) - + _ + ...\ . (33) 

'2^ ' 12e 720 £3 30240 ' ^ ' 



The analysis goes in a similar way. If one takes the third-order value Eq = 
0.95a;, one obtains |'(/'o(0)p ~ 0.9 In the e ^ oo limit one obtains the local 
duality relation 

{|^o(0)in^^ = y. (34) 

So, the "exact value" for the duality interval sq is 2uj. 

Summarizing, using the sum rule method one can determine the param- 
eters of the lowest resonance with « 10% accuracy. The main source of the 
error is the use of a rather rough model for the spectral density of the higher 
states. It should be noted that in QCD situation is better: normally only the 
first resonance is narrow, while higher states are broader and broader, so that 
approximating them by free quark functions does not lead to large errors. 



"^Plotting Eq for various values of is absolutely straightforward, and wc encourage the 
readers perform the fitting as an excercise. 
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2.2 Sum rules for form factors of the oscillator states. 
To calculate the form factors 



Fu{Q)^ rk{x)Mx)e''^^d''x (35) 



within the sum rule approach, one should construct a similar integral using 
the Green functions (llT 



G{x, l/^e|0,0) = X]1^^(a;2)1^fe(a;l)e-^^/^ (36) 

fc=0 

This gives the amplitude 

M(ei,e2,g2) ^ [ G{x, 1/zei | 0, ) G*( x, l/iea | 0, )e*0" 



e'^^ijl{x)^i{x)d^x 



g — -Efc/ei —El/ €2 



= J2 ^r(0)^fc(0)Ffc,(Q2)e-^^'=Ai-^'/^^ (37) 
fc,;=o 

One can write this amplitude in the (double) spectral representation: 

M(ei,e2,Q') = — / / e-^i/^i-^^/^V(si,S2,Q')dsids2 (38) 
Jo Jo 

where the spectral function is the sum of double delta-functions: 

oo 

p(si,S2,Q') E 4'!{0)M0)Fki{Q)S{si~Ek)S{s2-Ei). (39) 

k,l=0 

In the free case one has 

Gf'-^^ix, l/ie\0, 0) = !!^e-"-^^/2 ^ / J^^^kx ^-k- /2me (40) 

27r J (27r)^ 

for the Green function and 

M/'--(ei,e2,Q')= / ^e-'='/2'"'^-('=+'3)'/2'"^^ 

V27r/ £1 + 62 V 2m(ei + 62)/ 
10 



for the amplitude. The free spectral frinction pf^^^{si,S2,Q^) can be easily 
extracted: 



S2 - 



2m 



2m 



_ e{\V2msl + y/2ms^\ > \Q\ > |^/2mi^ - y^gma^l) 
The restriction on si , S2 is in fact the triangle relation 

|fcl| + |fc2| > IQI > llfcll -|fc2|| 



.(42) 



between the initial momentum fci (|fci| = \/2msi), the momentum transfer Q 

and the final momentum k2 = {ki + Q) (|fc2| = \/2ms2)- 
The oscillator Green function is also known: 



C'^ix, l/ie|0, 0) = 



27rsinh(a;/e) 
and using it, one can calculate the amplitude 



exp - 



mux^ cosh(a; / e) 
2sinh(a;/e) 



(43) 



M(ei,e2,Q') 



mo; 



27rsinh(a;/ei +a;/e2) 
Now, representing M as 



exp - 



sinh(a;/ei) sinh(ti;/e2) 
2mw sinh(a;/ei + w/ei) 



(44) 



M{ei,€2,Q^) 



TOW 



exp 



■ ^ g-(2fe+l)a;(l/£i + l/£2) 
fc=0 



\ 4mu) ) 



1 - 



(1 







/ X I (45) 
(1 - e-'^/«ie-'^/«2) \) ^ ' 

and expanding M in powers of e"*^/*^ = E\ and e"'^/^^ = E2 one can extract 
the form factors Fif.{Q^): 

M(e„e2,Q^) = ^E^i""^^r^expf-^) 



fc=0 



12 r 



1 



4mw 



1 



(l-£l)(l-£2) 



2 \4ma; / 



1 - S1S2 

identifying them as the factors in front of the S^S^ terms: 

^bo(<3^) = exp 
11 



+ ... 



(46) 



Amuj J 



exp 



Amuj 



■ 



Amuj 



J \ Amuj J 



(47) 



It is instructive to write down the amplitude M(ei, £2, Q^) for = 0: 



M(ei,e2,Q2 = 0) 



27rsinh(cLi/ei + ti;/e2) tt 



(48) 



fc=0 



It is easy to notice that 

• There are no terms with k ^ I: only the diagonal transitions are allowed in 
this limit 

• The numerical factor is just the same as in Eq. (|l^ ) for Af°''^(e). This means 
that i<fcfe(0) = 1 for all fc, just as required by charge conservation. 

Next step is to get the perturbative expansion in powers of w : 



M(ei,e2,Q2) 



27r(ei + £2) 



OJ 



1 1 



6 V^i £2 



X exp 



1 



2- (^ + i7)(i + f(^ + i)^ + ...) 



™eie2 
27r(ei + £2) 



ei £2 



1 + i 

6 \ei £2 



exp 



2m{ei + £2) 



To/ ^ 
360 €2 



6m£i£2(£l + £2) 



Q^cu^[ie,+e2r + 2e,e2] 



90to£3£3(£i + £2) 72m2£2e2(e^ + ^^)2 ■ 

= M/'--(ei, £2, Q^)\l + Y,{uj^f^k{^,,e2, Q^) 



The sum rule now reads 



(49) 



|^o(0)pFoo(g2)e-^«/^i-^°/^= + " 

= M/'-^^(£i,£2,g2)|i + o 



higher states'^ 
2 ^ 



O 



(50) 
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To proceed further, one should take some ansatz for the higher states 
contribution. The ansatz similar to that used earlier is to assume that 



^^higher states" = free states" 
outside the square (0, sq) (8) (0, sq). Then the sum rule has the form 

= ^ dsi ds2 p^'''"'{si,S2,Q^)e 

Jo Jo 



si/ei-S2/e2 



-o 



o 



In the ei, 62 ^ c» limit one obtains the local duality relation: 

|^o(0)pFoo(g') = — / / p^'-^^(5l,S2,Q2)dsidS2 
^ Jo Jo 

with p/'-^^(si, S2, Q2) given by Eq. @: 



^free 



{si,S2,Q^) 



d^kS si 



2m 



jk + Qf 
2m 



Using the two expressions given above, one can write 



1 



i2ny 



d^ke 



2m 



<so]0 



{k + Qf 
2m 



< So 



(51) 



(52) 



(53) 



(54) 



This integral can be easily calculated if one notices that it defines the overlap 
area of two circles having the same radius — 2mso , with their centers 
separated by the distance equal to Q: 



mso 



(k) - kVi - 6I(Q2 < Smso) , (55) 



where k — Q/^/SmsQ. Recall now that from the local duality relation for M{e) 
we had 

(56) 



{|^o(0)P}^^-"'° 



2tt 



and this gives the result 



Form-^ 

TT 



COS ^(k) — K\/l — 
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< 8mso) 



(57) 



having the correct normaUzation Fqq-'{Q'^) = 1. 

The local duality prediction calculated for the "exact value" sq = 2a; of 
the duality interval 

-i( Q 



FoO^ {Q'^)\so=2w — - 



cos 



VT6mw V ISmw 
is very close to the exact form factor 



6»(Q2 < 16mw) (58) 

F-*(Q^) = exp(-^). (59) 

There is a very simple physics behind the local duality. One should just 
compare the local duality prescription for the form factor 

'{k + Qf 



2m 



< So (60) 



and the exact formula 

J2 



Fm{Q^) = / MkWoik + Q) (61) 

to realize that the local duality, in fact, is equivalent to using a model form for 
the wave function in the momentum space: 

'^o''^^^^ K^^^'''^^'"'°^ ^^^^ 

or, taking sq = 2oj: 

iPl^^{k) = J— e{k'^ < 4mw), (63) 



to be compared with the exact wave function 



Mk) = 2J—exp(-^). (64) 

It is easy to establish that the model wave function correctly reproduces the 
simplest integral properties of the exact wave function: 



i>'''""=\k)d'^k = J i^^'^{k)cPk (65) 
j \i^^^'"=\k)f(fk = J \ip^^{k)f(fk. (66) 
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3 Condensates ant Operator Product Expeinsion 

3. 1 Correlators 

The basic object in the quantum mechanical example was the Green function 
G{xi, ti\x2, describing the propagation of a particle from one point to 
another. In fact, it was sufficient to analyze the probability amplitude for 
the partricle to remain at the same spatial point after some imaginary time 
elapsed: 

M(e) = G(0, 0|0. 1/ie) . (67) 

The question is what should bo its analog in the quantum field theory? 
Recall, that the motion in the external potential corresponds to a two-particle 
system with one of the particles being infinitely massive. If this particle moves, 
the other moves together with it. So, the Green function describes really a 
collective motion of two particles from one point to another. In quantum field 
theory this corresponds to the Green-like functions 

{Q\T{j{x)j{y)m 

normally referred to as correlators, with the currents j{x) being the products 
of fields at the same point, e.g. J^(x) = q(x)^^q{x) . Due to the translation 
invariance 

(0|T {j{x) j{y) )| 0) ^ (0|r - y) M )| 0) 

it is always possible and convenient to take one point at the origin. In QFT, it is 
also more convenient to consider correlators in the momentum representation: 

n(g) = j e'«-(0|r(j(x)j(0))|0)d4x. 

In the lowest order the correlator is given by the product of the relevant prop- 
agators, but there are also radiative corrections, which one can calculate in 
perturbative QCD, at least in the region where q is spacelike < and large 
enough. 

Next step is to relate n(g) to a sum over the physical states. This is 
performed via the dispersion relation 

niq^) = - f ^^^\ ds + ''subtractions" . 
71- Jo s - 9 

The subtraction terms appear if the spectral density is such that the dispersion 
integral diverges. If, say, p{s) ~ s", one needs (n + 1) subtractions: 
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ds p{s) 



A2 (s-A2)2 



(s - A2)»+i ^ 



+ {n(A2) + _ A2)n'(A2) + . . . + ,^'^" n(")(A^)l 

n! J 



(68) 



After collecting all the terms containing p(s), one observes that p{s) is divided 
by {s— q^){s~ X^)"^^^ ^ and the integral converges. Note, that all the subtraction 
terms are polynomials of a finite order in (from now on we use the notation 



3.2 Borel transformation 

Note, that in the oscillator studies, we dealt with the exponentially weighted 
sum: 



M{e) 



(69) 



We observed then that it has a nice property that for e values of an order of 
cj, both the contributions of the higher states were strongly damped and the 
perturbative expansion for A/(e) was converging fast enough. With the power 
weight implied by the dispersion relation, the suppression of the higher states 
contribution is not so effective. 

However, it is rather easy to transform the original dispersion integral into 
the exponentially weighted one. To this end, one should apply the so-called 
Borel transformation, formally defined by the formulaJ 



$(Af2) 



1 



lim 



(Q 



2\n 



n(Q2). 



(70) 



It is easy to check that applying the above procedure to l/(s + Q'^) one really 
obtains an exponential. Indeed, 

1 r; I 

(71) 



rfQ2 



Multiplying the result by (^^^ly one obtains 



1 



1 



1 



1 



(S + Q2)"+1 Q2/„(l_^s/Q2)n+l 

Finally, taking the rt — > cx3 limit one obtains 



AP (l + s/nM2) 



n+l ■ 



(72) 



M 



e "Z*^ . This means that 



$(M2) = B(g^ Af 2) n(Q 



ds 
M2 



p{s)e-'/^'^' + ''nothing'' (73) 
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where "nothing" means that all the subtraction terms have disappeared, since 
no polynomial can survive the infinite number of differentiations. Thus, the 
function <^{AP) is the QFT analog of M(e). 

In many cases one can perform the Borel transformation using the formulaiJ 

B{Q^ Af2){e-'^0'} = ^(1 - AM^). (74) 

As an illustration, let us reproduce the result obtained above : 

/•OC 

B / e-"(^+«')da 

/ e-""(5(l - aM^) da = —je~'/^^' . (75) 
Jo 

Other functions, the Borel transforms of which we will need, are those having 
a power behavior in \/Q^. In this case we have 

1 2 a""^ 

{Q^r Jo {n-l)l 

Thus, a power expansion for $(M^) over 1/A'P converges much faster than 
the original expansion of n(Q^) over l/Q^, since each term is now suppressed 
factorially. This looks precisely like the Borel improvement of a series, and 
that is why Shifman, Vainshtein and Zakharo'Jil referred to the transformation 
they introduced as "the Borel transformation". The representation for $(Af^) 
has the form of the Laplace transformation, and that is why some people call 
the sum rules obtained in this way "the Laplace transform sum rules". In 
fact, there is no contradiction: $(M^) is the Borel transform of n((5^) and the 
Laplace transform of p{s). 

To summarize, the Borel transformation produces two improvements: 

• higher state contributions in the spectral representation for $(Af^) are ex- 
ponentially suppressed, 

• its power series expansion in 1/M^ is factorially improved. 
3.3 Heavy-light system and the m ^ oo limit 

Using the basic formula, it is rather easy to establisl0 that in the limit when 
one of the particles becomes infinitely massive, the QFT function ^{AP) con- 
verts into its quantum-mechanical analog. Consider the lowest-order diagram. 



B- 
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The propagator of the heavy particle is 

1 



{q + — q'^ — rn? + 2qk + fc^ ' 

In the m — > oo limit all the bound state masses are equal to this mass plus 
some finite quantity. So, it makes sense to subtract this infinite constant. To 
this end one can represent the total momentum q as 

q = mv + I 

where v is the four-velocity. Then 

q^ -rn^ = 2m{vl) + 

The scalar product (yl) = E just has the meaning of the energy variable. There 
is only one more 0{m) term in the propagator: 

2{qk) = 2m{vk){l + 0{l/m)} . 

With this accuracy, one can write down the relevant Feynman integral as 

4 1 



I(m,v,E)= I —A^—-S{k) 



where is the propagator of the light particle. Now, applying the Borel 
transformation B{E e) which converts + (vk)) into ^e~'y^^^l'^ we get 

M{m,v,t) = B{E^t)I{m,v,E) = ^^ J e-^''^'^''S{k)(fk. 

It is easy to realize that the integral above is just the light particle propagator 
in the configuration space 

i.e., the Green function describing propagation in the imaginary time, since the 
velocity vector in the rest frame of the heavy particle, is oriented exactly 
in the time direction. 

Thus, we have demonstrated that the Borel transformed correlators arc 
just the QFT generalization of the imaginary time Green functions we studied 
in the quantum-mechanical examples. 
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3-4 Perturbative vs nonperturbative contributions in QCD 
Prom the technical point of view, the strategy now is 

• to calculate the correlator n((5^) in the deep spacelike region where one can 
rely on the asymptotic freedom of QCD, 

• apply the Borel transformation to the resulting expression to get a 

expansion for $(M^), 

• finally, using the dispersion representation for $(M^), to construct a sum 
rule. 

The question is, what are the contributions analogous to the co^/e^ correc- 
tions? 

The standard idea about the QCD potential V{r) is that it consists of a 
Coulomb-Iike part ^ Q:s(l/r^A^) calculable perturbatively at short distances 
and a nonperturbative ^ r long distance "confining" part: 

V{r) = y"Co";om6"(^) ^ V''°'^^r) 

The magnitude of the "Coulomb" effects is determined by the QCD run- 
ning coupling constant 

The standard modern estimate for the value of the QCD scale A is 

A ~ 200 MeV . 

As we will see, the parameters that really appear in QCD sum rules have 
dimension of {mass)"^ , and one should compare 

A^ ~ 0.04 GeV^ 

to the typical hadronic scales like p and nucleon masses (squared): 

ml w 0.6 GeV^ m% w 0.9 GeV^. 

To get such scales from A^ one needs unnaturally large numerical factors of 
an order of 25 or even 100. So, it is most unlikely that the O(as)-corrections 
can generate scales of an order of 1 GeV and , hence, the "Coulomb" part 
is practically irrelevant to the formation of the spectrum of light hadrons. 
Our problem now is to take into account the effects due to the long-range 
nonperturbative part of the QCD potential. 

Our study of the quantum-mechanical oscillator (which is the simplest ex- 
ample of a system with confinement) has demonstrated that the propagation 
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amplitude in the presence of the potential is different from the free Green func- 
tion. Furthermore, wc observed that the difference vanishes at short distances 
and that one can calculate exact Green function perturbatively, expanding in 
powers of the oscillator potential. In QCD, with two components of V{r) one 
can imagine a double series, both in v"'-^"'^''""^^" (r) and V'^°^^{r). There is no 
problem to arrange an expansion in y "Coulomb" fj.^ ^ _ ^j^jg ^jjj j^g^ stan- 
dard perturbation theory in the QCD coupling constant. On the other hand, 
the confining potential "l/'^°"^(r) in QCD is not even known. Moreover, the 
widespread belief is that it is determined by some essentially nonperturbative 
effects, i.e., by those which cannot be seen in perturbation theory. 
In such a situation, a possible way out is to proceed as follows: 

• to construct perturbation expansion in terms of quark and gluon propagators 
(this amounts to taking into account only the Coulombic part of the potential), 

• to postulate that quark and gluon propagators are modified by the long-range 
confinement part of the QCD potential; but the modification is soft in a sense 
that at short distances the difference between exact and perturbative (free- 
field) propagators vanishes. 

To formalize this statement, one can write the exact propagator 2?^^'*'=* (x) 
as a vacuum average of a T-product of fields in the exact vacuum Cl (corre- 
sponding to a confining potential) 

V^^-^^x) = {fl\T{^{xMO)m. 

According to the Wick theorem, one can write the T— product as the sum 

T{^{xMO)) = ^{xMO) + : vixMO) : 

^ V ' 

of the "pairing" and the "normal" product. The "pairing" is just the expecta- 
tion value of the T- product over the perturbative vacuum 

^{x)ip{0) = {0\T{<fi{x)ip{0))\0). 

^ V ' 

i.e., the perturbative propagator. By this definition, the normal product 
: ip{x)ip{0) : vanishes if averaged over the perturbative vacuum: 

(0| : fixMO) : |0) = 0. 

Thus, our assumption that 'D'^^'"=*(a;) ^ 'Dp^'^*{x) is equivalent to the statement 

(0| : <fi{xMO) : \fl) ^ 0, 

which is the starting point to calculating power corrections in QCD. To simpliiy 
the notation, in what follows we will write the Ihs of the above equation as 

(ifixMo)}. 
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3. 5 Condensates 

In the oscillator case the analog of {(p{x)ip{0)) is the difference between the 
exact Green function 

G°-(0,0|0, r/i)= "T, . 

ZTT smh(a;r) 

and the free Green function 

G/-^(0, 0|0, T/i) = ^ 

where r = 1/e is just the imaginary time variable. Note, that both oscillator 
and free Green function are singular for t ^ 0, but the difference 



m 
2^ 



C'^O, I 0, r/i) - G^''""{Q, | 0, r/i) 

(77) 



. 1 ^2 — wV^ - wV^ + . . . 

6 360 15120 



is regular at that point, and one can even expand the difference into the Taylor 
series in (wt)^. The expansion is possible essentially because the confining 
potential is regular (and even vanishing) at r = 0. 

In the QCD sum rule approach it is also assumed that the confinement 
effects are sufficiently soft to allow for the Taylor expansion of ( ip{x)^p{Q)) at 
a; = 0: 

(<^(0)y>(a;))=^ ^ --^^ {md,,...d,MO)) 

= ( + ^d^^) + ^^^^{^d^^d^,^) + ... (78) 

This is, in fact, the expansion of the nonlocal object {Lp{Q)Lp{x)) over the 
vacuum matrix elements of the local composite operators. Not all operators 
really contribute to the expansion above. Using the Lorentz covariance, it is 
trivial to find that 

( ipd^^p) = 0, 

(</'^/^i5a«2</5) = (79) 
etc., so that finally one arrives at the expansion in x^: 

(^(0M-)) = E(Tj ^^[(^^ • 



n=0 
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Thus, the modification of the propagator by the nonperturbative effects is now 
parametrized by the matrix elements of the composite operators like {ip{d'^)"ip). 
The examples in QCD are 

• (qq) referred to as the quark condensate, 

• (qD'^q), characterizing the average virtuality of the vacuum quarks, 

• the gluon condensate (GJ^j,G°^) , etc. Here = 9^ — igA^j^ is the covariant 
derivative and = {i/g)[D^,Di,] is the gluonic field strength. Note, that 
only gauge invariant composite operators should appear in QCD, i.e. each 
must be accompanied by the relevant A^. How this happens, will be discussed 
later. 

3. 6 Operator product expansion 

It is instructive to analyze the above expansion in the momentum representa- 
tion 

(^(o)^(x)) ^ j e^'^^v{k)d'k = (t) J i-k'rT^('^)d'k. 

(81) 

Note, that existence of {(p{0)(p{0)) implies that I?(fc), the nonperturbative part 
of the propagator, vanishes for large momenta faster than k^^, while exis- 
tence of ((^(0)9^93(0)) implies that 'D{k) < k~'^, etc. All matrix elements 
((/3(0)((9^)"93(0)) exist only if, for large momenta, the function I?(fc) vanishes 
faster than any power of say, like an exponential. Thus, the function 

'D(k) is concentrated in the region of small momenta. This corresponds to 
the basic assumption that at large momenta one should enjoy the asymptotic 
freedom. 

In other words, the exact propagator is represented as a sum of the pertur- 
bative 0{l/k^) part and the nonperturbative ("condensate") part that vanishes 
fast when the momentum k increases. Using such a decomposition, we obtain 
a modified diagram technique for the correlator functions. Consider, e.g., the 
lowest order diagram corresponding to the product of two propagators. After 
the decomposition it produces four diagrams (see Fig.||), where the dashed 
lines correspond to the nonperturbative part of the propagators. 

The first term (Fig.^) is just the ordinary purely perturbative diagram 
having ^ \og{q^) behavior for large q^: 

/f(^-^°^-^'^^)- 

Two other diagrams (Figs.p3,c), as we will see below, behave like 1/q^. This 
behavior is completely determined by that of the perturbative propagator: the 
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q-k q-n 

a) b) c) d) 



Figure 3: Operator product expansion for scalar one-loop diagram. 

whole momentum flows through the perturbative line where its propagation is 
not suppressed. In the last diagram (Fig.^), all possible momentum flows are 
strongly suppressed, e.g., if I?(fc), the nonperturbative part of the propagator, 
vanishes exponentially at large fc, so does this contribution also. 
Now, let us analyze the contribution of the second diagram: 

J<^^>(q2) = y'p(fc2)_^. (82) 

As we discussed earlier, the function T>{k^) is concentrated in the region of 
small momenta. So it makes sense to expand the perturbative propagator 
l/{q-kf in k: 

where 

{gfc}(")^{g^,...g^„}{A:^V..A;^"} 
is the product of two symmetric 

Q...lli...flj... _ Q...Hj...lli... 

and traceless 

tensors formed from the vectors q and fc, respectively. The symmetric-traceless 
tensors have a very useful property: 

F{k^){k^'K..k^'^}d^k^Q 
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for n ^ 0, because it is impossible to construct a nontrivial symmetric-traceless 
combination using only the metric tensors g^fij ■ To do this, one should have 
at least one vector, e.g., 

Thus, in our case 



= ly v{k^)e{e < q')d'k + J v{k')^eiq' < e) 

= —{(fiifi) + ^^quickly vanishing contribution" , (84) 
where the "quickly vanishing" term is given by the integral over large momenta: 

^^quickly vanishing contribution" = J 'D{k^)9{k^ > (f')d/'k 

where the function is strongly suppressed. 

Thus, it is very easy to calculate this contribution: the whole external 
momentum flows through the perturbative line, so one should just multiply 
the perturbative propagator by the condensate ((^93) factor. One can 

wonder, however, why the higher condensates (containing the 9^) did not ap- 
pear? To answer this question, it is useful to analyze this contribution in the 
configuration representation. The original diagram decomposition now looks 
like 

V{x, 0)P(0, x) = VP"^'' [x, 0)I?f (0, x) + VP''-'' (x, 0) ( ip{Q)ip{x)) 

+ VP^^' (0, x){ vixMO)) + { mvi^)) ( fi^MO)) , (85) 

and its second term can be written (modulo trivial factors like 27r) as 

{cp{OMx))/x\ 

Now, using the a;^-expansion for the ( (fi{0)(fi{x))- factor, we obtain 

l^{^iO)^(x)) = ^{^^) + ^^{ipiOMx))+0{x') + ... . (86) 

X X X Q 

In the momentum representation, the first term just corresponds to the ( jq^ 
contribution, while the second is proportional to <5^(g) and can be discarded 
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since tlie external momentum q is assumed to be large. Higher terms are pro- 
portional to the derivatives of and also can be ignored for large q. It 
is easy to realize, that the sum of all these terms just produces the "quickly 
vanishing contribution" discussed above. Only the terms singular for = 
(in the configuration representation) can produce a power-like behavior in the 
momentum representation. 

The expansion_we discussed above is a simplified illustration of the operator 
product expansioitj 

nj(a;)j(0)) = ^C„(x2)O„(0). (87) 

n=0 

It represents the T-product of two currents j{x),j{0) as a sum of the local 
(composite) operators of increasing dimension (in units of mass) . The coeffi- 
cient functions C„(x'^) characterizing the contribution of a particular operator 
are ordered in their singularity at = 0. From a simple dimensional analysis 
it follows that there should be a simple relation between singularity of the co- 
efficient function and the dimension of the relevant composite operator; e.g., 
if 

then 

On - (m')'°+", 

SO that after the Fourier transform one arrives at an expansion in powers of 
/i^/g^. There are, of course, more complicated contributions into the OPE, 
corresponding to more complicated original diagrams. 

Similar decomposition of propagators into perturbative and nonperturba- 
tive parts can be also performed in the QCD case. At this level of discussion, 
the only difference is that in QCD one has to deal with two types of fields: 
quarks q, q and gluons A. Consider the decomposition of a two-loop diagram 
(see Fig.|). 

• The first diagram (Fig.^a) corresponds to the purely perturbative contribu- 
tion (unity operator in the OPE language), 

• the second diagram (Fig.^) (in fact there are two such diagrams) corre- 
sponds to the {qAq)-type contribution, 

• the third (Fig.^) is evidently proportional to the (AA) vacuum average, 
accompanied by the one-loop coefficient function, 

• the fourth (Fig.^) diagram (there is also another such diagram) brings in 
the four-quark vacuum average (qqqq) . 
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a) b) c) d) 

Figure 4: Operator product expansion for QCD two-loop diagram. 

4 Operator Product Expansion: Gauge Invariance and Related 
Problems 

Looking at the operators discussed in the previous section, one can notice 
that some of them, e.g., {qd^q), {qA^q), {A^A^), etc., do not satisfy a very 
important requirement: they are not gauge invariant. And the product of 
gauge invariant operators in a gauge theory should have an expansion in gauge 
invariant operators Uke {qD^q) and (G^^C* 

4.1 Notations 

Let us introduce now the necessary notations. To begin with, 

Dn = d^- igA^ 

is the covariant derivative in the quark (fundamental) representation of the 
SU{3)c gauge group: 

with matrices (t")ab simply related to the Gell-Mann (3 x 3) A-matrices 

t" = A"/2. 

The indices A, B correspond to 3 possible quark colors:^, B = 1,2,3. The 
gluons have 8 colors described by the a-indices: a = 1,...8. The covariant 
derivative 

acting on the gluonic field contains the (8 x 8) cr-matrices 

A, = A^a" 
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characteristic for the gluonic (adjoint) representation of the SU{3) color group: 

(o'°)6c — —ifabc, 
where fate are the structure constants of the group: 

The gluonic field strength tensor 

can be related to the commutator of the covariant derivatives: 

G^u = -[D^,,D,] (88) 
G^.^-[D^,D,]. (89) 

4-.2 External field method 

A rather innocent-sounding statement that only gauge invariant operators 
should appear in the operator product expansion implies that there are nu- 
merous nontrivial correlations between the coefficient functions of non-gauge- 
invariant operators originating from completely different diagrams. For ex- 
ample, the operator (qd^q) appears when one expands the nonlocal operator 
(q{0)q{x)) in the simplest combination (see Fig.|^) 

(qiO)S{-x)qix)), 

while the operator (qA^q) necessary to form a gauge-invariant sum, results 
from a more complicated combination(see Fig.^) 

[mS{~y)ri^g)A^iy)Siy - x)qix)) d'y. 

In both cases, the perturbative propagators (in configuration space) are written 
explicitly, while the nonperturbative parts, in accordance with our previous 
discussions, are denoted by vacuum averages of the relevant normal products 
of fields. By a straightforward integration one can establish that 

S{-y)rS{y - x) d^y = -x^S{-x) , 
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a) b) c) d) 



Figure 5: Summation over gluons inserted into hard quark propagator. 



and, hence, the coefficient functions of the two operators are just those neces- 
sary to form the gauge invariant combination. 

However, this is definitely not the method one would like to use in more 
complicated situations, fn fact, it is much easier to analyze the sum of all the 
combinations differing only by the number of the gluons going into vacuum 



One can write it a£3'til {q{0)S{0,x; A)q{x)) , where S{0,x;A) has the 
meaning of the perturbative quark propagator in the (external) vacuum gluonic 
field A. It satisfies the Dirac equation 



Solving it by using perturbation theory in A, one gets back the original dia- 
grammatic expansion. 

The function S{x, y; A) is not gauge invariant. From the good old electro- 
dynamics it is known that propagation of a charged particle in the external 
field results in a phase factor given by the exponential of a line integral. In the 
QCD case it looks like 



Jc 

where P means that the color matrices should be ordered along the path C 
connecting the points x and y. Furthermore, it is known that the combination 



is gauge invariant (though path-dependent). This gives the idea how one can 
construct the operator product expansion in an explicitly gauge invariant way. 
Let us try the ansatz 



(see Fig.|). 



^ '^d^^ix)) - m S{x, y; A) = -S^{x - y). 



(90) 




{q{y)P{y,x-A\C)q{x)) 



S{x,y-A)^E{x,y-A)S{x,y-A) 



(91) 
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where E{x, y; A) is the P-exponential for the straight-Hne path connecting the 
points X and y: 



E{x, y; A) = Pexp 



z=y+t{x-y) 



It is straightforward to derive that the original Dirac equation is satisfied 
only if the function S{x, y; A) is a solution to the modified Dirac equation 
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igA^{x;y) ) - m 



S{x,y;A)=-5\x-y) 



that differs from the original one only by the change 



Al{x)^Al{x;y) 



= {x''-y") f 
Jo 



GUz)E''^z,y)tdt 



(92) 



(93) 



Here, z = y + t{x — y) and E \s & straight-line-ordered exponential with the 
j4-fields multiplied by the cr-matrices of the gluonic representation. It appears 
after one commutes the £^-exponentials through a r-matrix: 



{T'')ABEBc{z,y) = EAB{z,y){T^)BcE''\z,y). 

The commutation rule is based on the well-known formula 

1 



and the relation 



e^Be-^ = B + [A,B] + - [A, [A, B]]+... 



[To, Tb] = —{(Ja)bcTc- 



(94) 



(95) 



(96) 



Surprizingly enough, the presence of the £^-factor only simplifies the situ- 
ation, because the derivatives in the Taylor expansion of the G£^-combination 
are just the covariant ones: 



oo 

Gl^{z)&^{z,y) = ^ - G^^;^,...^„(j/)(z - yT' ...{z~ yY" . 



(97) 



n=0 



As a result, the new field A'^{x; y) can be expressed in terms of the gluon field 
stress tensor G^,/ and its covariant derivatives (taken, of course, in the gluonic 
representation): 



A''^{x;y) = Y, 



1 



n=0 



{n + 2)n\ 



{x - yfix - yY- ...{x- yY- G»^^^^...^Jy). (98) 
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This means that the function S{x, y; A) depends on the giuonic field only 
through the gluon field stress tensor G^jy and its covariant derivatives. Solving 
the equation for S{x, y; A) perturbatively, one obtains the (;^-expansion that 
has the same structure as the original (7A-expansion. A very important dif- 
ference is that the resulting operators have an explicitly gauge invariant form. 
For example, taking the lowest (zero) order term one obtains the contribution 

{q{y)S{y-x)Eiy,x;A)q{x)) (99) 

containing the gauge-invariant operator 

oc ^ 

{q{y)E{y, X- A)q{x)) = ^ " V^^^^ ' ' ' " 2/)m„ (9^^^ • ■ • ^^"9> \y (100) 

n=0 

We used the direct straight-line path from a; to j/ in our ansatz. In general, 
it is more convenient to use the path formed by two straight lines: from x to 
some fixed point zq (usually taken at the origin, zq = 0) and then from zq to 
y. The only change then will be that quark and gluon fields q{C)^ ^(0 '^iU be 
finally accompanied by the i?(^, zq\ A) exponentials, so that a covariant expan- 
sion will be obtained if one expands at the fixed point zq. Using translation in- 
variance, one can always substitute matrix elements {q{zo) . . . G{zq) . . . g(zo)), 
resulting from such a procedure, by (g(0) . . . G(0) . . . q{0)). Less trivial is the 
disappearance of the ZQ-dependence from the coefficient functions. One can 
enjoy it only if the calculations had no errors, and after summation over all 
diagrams (of the modified expansion) producing the same operator. The latter 
observation suggests that zq works as a gauge parameter. Indeed, it is easy to 
notice that the relevant P-exponentials 

E{x, Zq; a) = Pcxp[ig{x^' - z^) / A^(z) dt]\2=,^„+t(^^_,g) (101) 

Jo 

disappear (are equal to 1) if one imposes the gauge conditioiEl'El'00'0 

{x^-z^)A^{x)^0. (102) 

Furthermore, the "calligraphic" field A in this (Fock-Schwinger or fixed point) 
gauge coincides with the ordinary one A = A{{G}), i.e., in this gauge not only 
G^i/ can be expressed in terms of A^, but also the vector-potential A^ can be 
represented as a (nonlocal) functional of the field strength tensor G^j/. This 
means that theuEock-Schwinger gauge is a physical gauge. Another example is 
the axial gaugft^ 

n^A''{x)=0 (103) 
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with being a fixed vector. In this gauge one can also expresS in terms 
of G^^: 

A^{x)=n'' I G^f,{x + tn) dt. (104) 

^0 

Imposing the Fock-Schwinger gauge condition, one can forget about the 
exponentials. One should remember, however, that the "ordinary" derivatives 
in this gauge should be treated as the covariant ones, e.g., they do not com- 
mute, etc. To illustrate the typical tricks and conventions, let us consider the 
simplest gauge invariant nonlocal condensate 

{qiO)EiO,x;A)q{x))\^^A)=o = imii^)) 
oo 2 

= E -^/^i ■ ■ ■ {qD^' ■ ■ ■ D'^-q) = im) + ^{qD^) + ■■■■ (105) 



The ratio 

^ (106) 

m) 

can be interpreted as the average virtuality of the vacuum quarks. It is possible, 
however, to represent the (qD^q) operator in a different form: 

{qD\) = {qD^^D^q) - {qD^^ D'' q)g^, = (qppq) - {qD^^ D'' a^^q) 
= -mliqq) - ^(q[D^^ , D"-] a^^q) = -ml{qq) + ^iqigG^" a ^,q) (107) 

where we have used the identity 

7m 7^ - 7.^ 7m 
Qt^v = l^ilu ^ = It^ly - cr^^, 

the equation of motion 

$)q{x) = -imqq{x) (108) 
and the definition of the field strength tensor 

[D^,I),]=-igGf,,. (109) 

Thus, the "average virtuality" of the vacuum quarks (qD^q) is directly related 
to the "average vacuum gluonic field strength" {qigG^'^a^uq) . In many papers 
one can find the notation 

{qigG^'''a^,,q)=ml{qq). (110) 
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Using it, one can write the relation 

X' = ^-ml. (Ill) 
For light quarks, the term can be neglected. 
4-3 Equations of motion and the current conservation 

The equation of motion Ipq{x) = —imqq{x) brings in correlations between 
matrix elements of different operators. These correlations, in fact, are very 
important to preserve the symmetry properties required by the structure of 
the currents j{x) in the correlators. For example, one might expect that a 
correlator of the electromagnetic vector currents 

^,..{p) = / (T j^(0)j,(x))e^f- d^x (112) 

should be transverse: 

n^.(p) = (P^P. - 5^.p') n(p2) (113) 

due to the current conservation condition — 0. Indeed, for the pertur- 

bative one-loop contribution, using the explicit form of the propagator in the 
configuration space 



S{x) oc -2i 
one finds that the integrand is 



x"^ 



1 m 
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+ ^+0(m3) (114) 



[g^^x'' - 2x^x^] - [ff^,a;2 - 4x^x,] + . . . . (115) 

It is easy to see that both terms satisfy the transversality condition 

and this means that 

nr;*(p) = (p^p^ - 5^.p')n^'-*(p^) 

at least up to the third order in m. One can, in fact, perform the calculations 
in the momentum space directly, to see that this property of n^^^*"* (p) is valid 
to all orders in m. 
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{q{0)lf^S{x)j.q{x)) cx (9(0)7^ ( -2ij^ + — ] j,q{x)) + 0{m'') 



Now, let us consider the simplest nonperturbative contribution. In the 
momentum representation it looks like 

{11 ^i— 2^vq) ^{m) ■ (116) 

However, it seems that there is no way to get the Pf^p,^ term, since there is 
only one term in the numerator. To solve this puzzle, let us consider this 
contribution in the configuration representation. Up to 0(m2)-terms one has 

= -j^imi^^M^)) + ^imif^iMx)) + o{m^). (117) 

Note, that vacuum averages of local operators with an odd number of indices 
vanish, and only those with an even number survive. Incorporating the Lorentz 
invariance, we obtain 

(9(0)7/^71/9(0)) = 9f,u{qq), 

(9(0)7/*7a7i'-D/3g(0)) = {9t,a9i^a + QvaQua - 9iJ.v9au){qiaDl3q) , 

{qj,D0q) = ^{qPq). (118) 
Finally, incorporating the equation of motion 

(#(?) = -im{qq) (119) 
one obtains the transverse result 

{Q{OhiJ.Six)'y„q{x)) oc ■^^^{qq)[9nux'^ + 2x^x„]. (120) 

It is instructive to analyze also the (GG) contribution. First, using the 

exponential representation, it is easy to see that for the closed loop the ex- 
ponentials corresponding to the upper and the lower lines cancel each other, 
and only the 0{G) terms remain. This general result can be illustrated on 
the example of one-loop diagrams with two vacuum gluons. To study, how the 
cancellation proceeds in this case, one should expand the exponential up to 
0{A^) terms: 

S{x,A)=S{x)ll + ig{xA)+^-^{xAf 



{i9) 



2 



S{-x, A) = S{-x) { 1 - ig{xA) + ^{xAy + ...}. (121) 



■I 
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Now, combining the contributions of the three relevant diagrams, one obtains 



{igfl^l(xAr + ^{xAf - {xAf^ = 0, (122) 

i.e., the A— dependence disappeared. To calculate the 0{G)- terms, it is most 
convenient to use the Fock-Schwinger gauge (xA) = 0, in which 

M^) = G.f^itx) tdt = (n + 2)n! 

= ^x^G,^{Q) + .... (123) 



Using this representation, one can construct the Feynman rules for the vertices 
corresponding to G^^, G^v;^, etc. For example, in the momentum represen- 
tation 

2 ^""^ 2i''^dK' 



T a, ^ tP^.j^^ (124) 



and the insertion of such a vertex into the quark propagator gives 

= 1.1.1 (125) 

The result looks like an insertion of two A- fields, A,^ and A^, antisymmetrized 
with respect to the interchange of /x and u. There is an important thing 
to remember while constructing the diagrammatic technique for the G-type 
vertices: in the momentum representation they contain derivatives acting onto 
all the propagators between the G-vertex and the Fock-Schwinger zero point. 
The G-inserted propagator looks simpler in the configuration space: 

^(^. G) = -^J^ - x^e^p,M^^^p^, + 0{zo){GG) + ... (126) 

where zq is the Fock-Schwinger fixed point. If one takes zq = 0, then the 
0({GG)) term disappears. This means that in the gauge {xA) = only the 
diagram with the G-insertions into diff'erent lines survives. Its contribution in 
the configuration space is easily calculated by just multiplying the two S{x, G) 
propagators. Taking the trace over the 7-matrices, one obtains the result 

2i 

2\27^(5m^^^ + '^x^,Xu){g^G'^) (127) 



3(167r2)2(x2) 
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proportional to the same transverse structure as the contribution due to the 
mq{qq) operator. This is because both operators have the same dimension, 
equal to 4 in mass units. 

The QCD sum rules for the "standard" mesons, like p, tt, K, etc., are based 
on the OPE including the operators of dimension 6, at least. One of these 
operators can be easily obtained from the Taylor expansion of the nonlocal 
quark condensate containing a 7-matrix 

{miMx)) = {qiaDaq)x" + ^{qj,Da,Dc.,Dasq)x"'x"-x"' +... . (128) 

We recall that the terms, containing an even number of derivatives, vanish since 
the total number of indices is odd in these cases. Incorporating the Lorentz 
invariance, we obtain 

+ {qb''^bo) + {qb^b^p)). (129) 

Next step is to use the translation invariance {d^q{0) . . . q{0)) = , the equa- 
tion of motion .^q' = —irUgq and the commutation relation [Da,Dij] = igGa^ 
to get 

mhMx)) = ^^(^{qq) + ^{qb^q)^+^x'^x^q^,[b^,Gnq). (130) 

The commutator [£)a,G'*"] is equal to the covariant derivative of the gluonic 
field 

[^„,G''"] = G"'«;„ . (131) 
The final step is to use the gluonic equation of motion 

G"'«;a=i^"= J2 iqirraqi)T'^ (132) 

i—u^d.s 

analogous to the Maxwell equation. Thus, due to the equations of motion, the 
apparently two-quark operator produces a four-quark term: 

rnhMx)) = ^^((9-g) + ^imb'q}) + ^x<^a;2(MT"5). (133) 
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5 QCD Sum Rules for p-meson 



5.1 Dispersion relation 

As pointed out by SVZ (whose analysifl we closely follow in this section) we 
closely to determine the basic parameters characterizing the p-meson, e.g., its 
mass TOp and its e~'"e~ coupling constant gp, one should consider the correlator 
of two vector currents 

]^^^ ^\{ui^u-d^^d) (134) 
with the p-meson quantum numbers: 

(135) 

Via the dispersion relation 

nf") {q' = -Q') = nf") (0) r ds (136) 

127r^ 7o s(s + Q-) 

this correlator is related to the total cross section of the e~'"e~-annihilation to 
hadrons with the isospin 1 = 1, measured in units of the e+e^ fJ-^ fJ-^ cross 
section : 

n/=i a {e+ e- hadrons,! = 1) 
(T(e+e 11+ fi ) 

Indeed, the jjf'' coincides with the isovcctor part of the electromagnetic cur- 
rent and is responsible, therefore, for the production of the hadrons with 1=1 
in the e+e~-annihilation. The correlator nji^^(q) contains, of course, only 
the hadronic part of the relevant factors, and that explains the division by 



5.2 Operator product expansion 

On the other hand, one can construct the operator product expansion for the 
correlator n|i^^((7) calculating the relevant Feynman diagrams. To begin with, 
one should calculate the perturbative diagrams shown in FigJ^ to extract the 
perturbative spectral density 

pP-*(^)^3^(^^^ {l + ^} (138) 

which is essentially the isovector (1=1 isospin projection) quark charge 
squared (|)^ multiplied by 2 {u- and d- quarks) and multiplied by the number 
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a) 



b) 



c) 



d) 



Figure 6: Perturbative contribution up to two loops. 




Figure 7: Nonpcrturbativc contributions up to dimension six. 



of colors (3). The simplest nonperturbative contribution is due to the quark 
condensate (see Fig.^). It is proportional to the current quark masses: 



|m„MM + niiiddlO) 



with niu ^ AMeV, ~ IMeV . Hence, this contribution is rather small 
numerically. The contribution 

CG 2 (0|gG|0) 

due to the gluon condensate is given by three diagrams shown in Fig jTf^-d. As 
mentioned at the end of the previous section, only the diagram |^ contributes 
in the Fock-Schwinger gauge Xfj_A^^{x) = 0. It is rather straightforward to 
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calculate the diagrams |^e,f proportional to 

However, the contributions of the same four-quark type can be obtained also 
from the diagrams 0g,h which originally are bilinear in quark fields. Notice, 
however, the numbers 1 and 3 attached to the relevant lines going into the 
vacuum. They indicate how many covariant derivatives are acting on the 
relevant "vacuum" field. Thus, the matrix elements are {q'^iiDnGapq) and 
{q^^Du^Du-iDy^q) ■ As discussed at the end of Section 4, after vacuum aver- 
aging and commutations only the {q'y^D,jG f^^q) component survives and, by 
equations of motion it is reduced to (0|g7^r°j°(7|0) with 

jt= E (0|ft7MT>.|0). 



5.3 Basic QCD sum rule for the p-meson channel 

After the Borel transformation that modifies the dispersion integral 

(139) 



(140) 



s + Q^ J ' M"^ 



and the power corrections 

1 



(Q2)n (n-l)!(M2)'^ 
one obtains the following QCD sum rul(0: 



Jo 2 



as(M) 4:7:^,, , , T«^ 

TT M 

-^{as{u^^T-u + d^^T^d) J2 (019.7^^'^ 9.10)^141) 

i—u,d,s 



It is the starting point for further analysis. 
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5-4 Fixing the condensates 



Next step is to specify the numerical values of the matrix elements on the r.h.s. 
of this sum rule (the same values will be used in all other sum rules). 

• The first matrix element is in fact fixed by the current algebra: 

(m„uu) + (mjd) = = _i.7 . lO-^GeV^ (142) 

• The value of the gluonic condensate was extracted from the analysis of the 
QCD charmonium sum rulefl 

{^GG) = 1.2 • lO^^GeV^ (143) 



• The four-quark operators are first approximated by the value dictated by the 
vacuum dominance (or factorization) hypothesis: 

{'4'{A}'4'{B}'4'{C}^{D}) = {lP{A}1p{B}){'il>{C}i^{D}) " (V-jA} V'{D}) ('/'{C}V'{B}) 

(144) 

where the subscripts A, B, C, D include spin, color and flavour. Using the 
orthogonality relation 

{^{A}^{B})^^^m (145) 
one can reduce then any four quark operator to the (qq)'^ form. In our case 

vacuum dominance — 

32 

= yaaiw)^ ~ 6.5 • 10"'' GeV^ (146) 

((^Tck'^ ^) ^ ^ {Qi'ya^ Qi) {vacuum dominance — 
i—u.d,s 

-6.5- 10-4 GeV^ (147) 
with the numerical values following from the estimate 

asiqqf ~ -1-8 • lO^^'GeV^ (148) 

which can be justified using the above relation between the quark condensate, 
quark masses (m„ + = 11 MeV, experimental values of the pion mass and 
pion decay constant under the assumption that the relevant normaliza- 
tion point corresponds to = 0.7. In fact, this estimate amounts to fixing 
one of the most important input parameters of the QCD sum rule method 
and the quoted value is supported by successful predictions for many different 
resonances. 
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5.5 Qualitative analysis of the QCD sum rule 

Using the vacuum dominance hypothesis, we can rewrite the sum rule in a 
more compact form: 



1 



Af4 



TT .a. 



GG) 



3M4 ' TT ' M6 
Substituting the numerical values of the condensate parameters gives 



(149) 



/ e-'^''" R'=\s)ds = U'l' 
Jo 2 



1+— — ^+0.1 



0.6GeV^ 



-0.14 



0.6GeV' 



M2 
(150) 

Our experience with the sum rules for the quantum-mechanical oscillator (see 
Section 2) suggests that to calculate the p-meson mass we should incorporate 
also the "daughter" sum rule resulting from the original one after the differen- 
tiation with respect to l/Af^: 



as(M) ,0.6GeV^, .O.eCeV^^. 



TT 



Af2 



(151) 

Note now that power corrections in the sum rules above are small for the values 
of the Borel parameter M as low as AP = = 0.6 GeV^: 

/>CXD o 

/ e-'''<R^=\s)ds = -mill + 0.1 + 0.1 - 0.14], 
Jo 2 

/>oo o 

/ e-'^"''sR^=\s)ds = ^m^[l + 0.1 - 0.1 + 0.28]. 
Jq 2 



(152) 
(153) 



At such a low value of ]VP one should expect that the integral over the physical 
cross section is dominated by the lowest resonance, i.e., by the p-meson. Its 
contribution to i?^^^(s), in the limit of the vanishing width, can be written as 



K i') = ^S{s-ml). 

9p 



(154) 



These observations allow one to make a rough estimate of the predictions 
following from the QCD sum rule. To this end we, following SVZJl 
• neglect, for = nip, the power corrections in the r.h.s. of the sum rules, 
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• neglect, for A-P = m^, the higher state contributions in the l.h.s. This gives 
the following relation 

— ^e-i«|m^. (155) 

Simplifying it results in a famous SV^l expression for the gp constant in terms 
of fundamental constants tt and e: 

^ - ^ « 2.3, (156) 
in e 

to be compared with the experimental value 2.36 ± 0.18. 



5. 6 Fitting the spectrum parameters 

The standard procedure is to extract the parameters of the hadronic spec- 
trum requiring the best agreement between the two sides of the sum rule. Of 
course, having only a few first terms of the operator product expansion one 
should not hope to reproduce all the details of the resonance structure in a 
particular channel. Only the gross features of the spectrum can be extracted 
from the QCD sum rules. So, let us take the ansatz "first narrow resonance 
plus continuum" which has been discussed in the oscillator case. In fact, one 
should expect that in QCD this ansatz should work better, since the higher 
hadronic resonances are very broad, and the total cross section of their produc- 
tion is known to be well reproduced by the parton model (i.e., perturbative) 
calculation. Specifically, we take 

K^') = + + v) 0{s ~ so). (157) 

The parameters to fit are the coupling constant g^, the p- meson mass and 
the effective continuum threshold sq. 

The strategy will be just the same as used for the oscillator states: 

• transfer the continuum contribution to the right hand side, 

• write the original sum rule in the form 

y2g-m^/Af^ = ^Ap (l - e"""/*^') + ''powcr corrections'' , (158) 

where, for brevity, we introduced the notation — ^'^'^ ™p ^ 

• write the daughter sum rule in the form 

flmle-^'"^" = ^M^ (l - (1 + ^)e-''«/^^')-f "power corrections", (159) 
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• divide the second sum rule by the first to get the expression for as a 
function of sq and Af^ (for fixed values of the condensates), 

• look for the value of the effective threshold sq that produces the most stable 
curve for (recall that with all power corrections taken into account, and 
the exact ansatz for the higher states, one should obtain an A/^-independent 
result for m^). 

In this way one would obtain]^ sq ~ 1-5 GeV^ and rn^ w 0.58 GeV^. Next 
step is to extract from the first sum rule, using these values of sq and m^. 
The value sq ~ 1.5 GeV^ provides the most stable curve for g^, corresponding 
to Ri 2.17. The result is rather sensitive to the m^-value. If one takes the 
experimental value « 0.6 GeV^, then the most stable curve still is that for 
So ~ 1.5GeV^, but the t^^-value is a little higher g^ « 2.40 and closer to the 
experimental one g^ « 2.36 ± 0.18. Anyway, it should be emphasized that all 
the 5^-values extracted from the QCD sum rule differ from the experimental 
value by less than 10%, the claimed precision of the method. 

6 QCD Sum Rules for Pion 

6.1 Axial vector current 

Vector currents are a special case, because the relevant correlators can be re- 
lated to directly measurable cross sections. From purely theoretical point of 
view, other currents such as scalar, pseudoscalar, axial vector etc., are not 
worse, but experimental information in these cases is normally limited to the 
masses of the lowest states. Still, the axial vector current is in a better situa- 
tion. First, its projection onto the pion state is just the pion decay constant 

(0|J757m«I^) - (160) 

Second, measuring the decays of the heavy lepton r Vr+X, one can measure 
the coupling of a hadronic state X to the axial current. Our goal now is to 
show how one can use the QCD sum rules to study the hadronic spectrum 
in the axial vector channel. Again, our presentation here closely follows the 
classic SVZ analysiil. 

Since the axial vector current a^(x) — ^757^1* is not conserved, the relevant 

''With all the numbers explicitly given, the readers are encouraged to perform the fitting as 
an excercise. 
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correlator Il\^J^'{q) is a sum of two independent functions: 

K'u\Q) = i J e*«^(0|T(a+(x)a;(0))|0)d*a: = -5^.ni(Q2) + g^g,n2(Q2), 

(161) 

or, in terms of the transverse and longitudinal components 

n^rH?) = - [9,^ - ^} ni(Q2) - ^{ni(Q2) + q'mq')}. (I62) 

6.2 Massless pion 

The common bchcf is that the pion mass vanishes in the hmit of exact chiral 
symmetry so that the pion is a Goldstone particle. It can be shown, that 
the operator product expansion in this case really requires the presence of a 
massless particle in the limit rUq —^0. In this limit the imaginary part of 
the longitudinal function Ili{Q^) + Q^Il2{Q^) = ^\\ vanishes, and XIn is a 
polynomial in Q^. Now, let us switch on a small quark mass and keep terms 
linear in this mass. The statement is that the function 'n.i{Q^) + Q^Il2{Q^) is 
exactly calculable in this approximation; 

TT fr,2\ , ^2tt fr,2\ {mu+md}{uu + dd) 2s 

UiiQ^) + Q^U2{Q^) = ^ '- + 0{m^g). (163) 

There are no terms linear in nig of higher orders in l/Q^. Let us show that 
this equation really implies the existence of a nearly massless pion. First we 
should express the Il-functions in terms of the hadronic contributions, using 
the dispersion relations: 

ni(Q2) = - / -^ds + ''subtr." (164) 

1 1"°^ Tm TT 
n2(g2) = - / -^ds + ''subtr." (165) 

MQ') + Q^MQ') = - r i"^ni-5imn2 ^^ ^ „ ^^gg^ 

where Imlli contains contributions only from pseudovcctor (spin-1) states and 
their contribution should be transverse, i.e., proportional to ((?^i^ — q^qv I Q^)'- 

Imni = ^7r/l<5(s-M2), (167) 

A 
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while Inin2 contains contributions both from pseudovector (spin-1) states and 
from the pscudoscalar (spin-0) ones and their contribution into 11^"^^ (q) should 
be longitudinal i.e., proportional to q^qi,: 

slmHa = J2 " ^1) + Yl -^flmp^s - Ml). (168) 

A P 

Hence, only pseudoscalar states contribute to the longitudinal component of 
nltt^^ For small quark masses we have 

{mu + md){uu + dd) . ^, 2\ _ /p"ip _ 



Q2 ■ V ^Q2+^2 



AU 0((5~^) terms should be 0{m'j^). This means that 

flm% = 0{ml) (170) 

for all states. Hence, the particles with masses that remain constant as niq — > 0, 
decouple in the chiral limit: 

ml = 0(m°), (171) 
/P - O(m^). (172) 

The last relation also means that they do not contribute to the term linear in 
nig. The only way to get this term is to assume that there exists a state for 
which the 0(mg)-term is brought by the mp-factor: 

m% = 0(mj), (173) 
fp = 0{m°). (174) 



Such a state is naturally identified with the pion. In addition, we get the 
well-known relation 

f^ml = -{mu + md){uu + dd) (175) 

demonstrating explicitly the vanishing of the pion mass in the chiral limit. 

To derive the starting statement, one should consider qfj,quIl''i^J'\q). By 
equations of motion it is related to the correlator of the pseudoscalar densities: 

= -i j e*«'^(0|T {{u{x)^r,d{x)) (J(0)75u(0)) ) | 0) d^x + const. (176) 
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The constant on the right-hand side is due to contact terms which normaUy 
arise if one differentiates a T-product. This constant corresponds to the 
term in the combination ni(Q^) + Q^^2{Q^) ■ Direct calculation of this term 
gives the result stated in the beginning of this subsection. 

Thus, the operator product expansion imphes that — > in the ruq — > 
hmit, and , hence, the pion mass should not be fitted from the sum rule. 
Within our accuracy one can set — from the start. 



6.3 QCD sum rule for the axial channel 

The operator product expansion for the correlator of two axial currents is com- 
pletely analogous to that of two vector currents, the only difference being the 
presence of the extra 75-matrices in the vertices corresponding to the currents. 
For massless quarks, however, this produces no changes both in perturbative 
and 0{{GG)) contributions. However, the four-quark terms shown in Fig.|^,f 
change sign, because one of the 7-matrix type terms is substituted by the 
condensate factor {qq) which has the structure of 1 with respect to the Dirac 
indices. 

After the Borel transformation the sum rule looks as follows: 



Att'^ 

q—u,d,s 



(177) 



In the axial channel we have two states that are the lowest among other 
states with the same quantum numbers: a pseudoscalar (the pion) and an 
axial vector (Ai-meson) ones. Thus, one can try the ansatz with two low-lying 
states ( "pion+y4i+continuum" ): 

Imn2(s) = nfjis) + 7r/l/(s _ m^J + i- (l + ^) e(s > 4^^^ (178) 

One can also treat the Ai as a part of the continuum. This, more crude 
approximation, corresponds to the simplified ansatz 

ImU^is) = nfJis) + ^ (1 + ^) e{s > 4^^^). (179) 

The constants are defined as follows: 

(0|d757pw|7r) = i/^P^, {0\d-f5l,,u\Ai)=ifA,e^, (180) 
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where is the Ai polarization vector. 



6.4 Fitting the sum rule 

Substituting the numerical values of the condensates into the basic sum rule 



gives 



'e-^/*^'lmn2(s)ds= ^ 



as{M) 



+0.1 



0.6GeV^ 



0.22 



0.6GeV 

M2 



TT 
2x 3, 



(181) 



The "daughter" sum rule is 



poo 

/ e-^/*^'lmn2(s)sds 

JO 



M2 

47r 



as{M) 



M2 



TT 
2x 3. 



M2 



(182) 



Note, that in the rn^ = approximation the pion does not contribute to the 
daughter sum rule, and the spectrum structure for sIm 112(5) = Imni(s) is 
more like in the vector channel: first non-zero mass resonance + continuum. 
This is reflected in the structure of the relevant sum rule: power corrections 
are negative. They are considerably larger than those in the p-channel, and 
one should expect that the Ai-mass (squared) is much larger than m^. 

The fitting procedure ^ is most straightforward in the case of the simplified 
ansatz. One should transfer the continuum contribution to the right hand side 
to get /2 as a function of s}^"^ and the Borel parameter M"^: 



47r2/2 = Af2 



1 - e' 



' /M' 



3M- 



(— GG) 

TT 



7047r^ 



.{my 



M4 

(183) 

t2 „„j A„2f2 _ 0.65 GeV^ 



The most stable curve corresponds to Sq ~ 0.7 GeV and 'Att 
the experimental value being 47r2/2 = 0.67 GeV^. 

One can also try the ansatz in which Ai is treated as a (narrow) resonance. 
The strategy is just the same as in the p-case. To this end, one should 

• differentiate the daughter (second) sum rule with respect to 1/M2 to get the 
third sum rule, 

• transfer the continuum contribution to the right hand side in all sum rules. 



Again, we encourage the readers to obtain the curves themselves. 
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• extract m^^ and Sq ^ from the ratio of the third sum rule to the second, 

• using the values m\ w 1.6 GeV^ and So'^^' ~ 2.4 GeV^ obtained from the 
previous fitting, find f\_^ from the second sum rule, 

• using the value 47r^/^^ 1.95 GeV^ obtained from the previous fitting, find 

from the first sum rule. The value « 0.69 GeV^ obtained in this way 
is very close to the experimental one 

j2^exp ^ 0.67 GeV^ The Ai mass value 
« 1.6 GeV^ extracted from the QCD sum rule also is in a very good 
agreement with experiment. 

7 QCD Sum Rules for Pion Form Factor 

7. 1 Three-point function 

To analyze the pion form factor^azE consider the three-point function, i.e., the 
correlator of the three current^iJO: 

T^M'P2) - ^ J e-^P^^+^P^y{{jpiy)J>^iO)j + {x)})d^xd\ (184) 

where 

Jf = - ^dj^d (185) 

is the electromagnetic current and 

ja = dislaU (186) 

is the axial current used in the preceding section to study the static properties 
of the pion. The pion contribution 

(Olj^y) \P2) {P2 1 J'-iO) \pi){pi \j+ {x)\0) (187) 

into this three-point correlator can be obtained by inserting the pion state(s) 
between the currents. The form factor F^{q'^) appears in the middle matrix 

{p2\J^mPi)^ip't+pt^)FAq'), (188) 

while the other two are projections of the axial current onto the initial and 
final pion states 

(O|J0(O)|P2) =*A(P2)/3 (189) 
(pi|j+(0)|0) =*A(pi)„. (190) 

The three-point function Tl^p{pi,P2) is the sum of different structures each 
characterized by the relevant invariant amplitude ti{pl,P2, q^). The first idea is 
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to study the pion form factor analyzing the invariant amphtude corresponding 
to the structure (p2)/3(pi)a(Pi +^2)- However, there are other structures 
((Pi)/3(p2)a, (pi)/3(pi)a, (p2)/3(p2)a that coiucidc with (p2)/3(pi)a in the hmit 
q ^ 0. This comphcation disappears if all the basic structures are expanded 
in P = (pi +p2)/2 and q = p2 — Pi- Then the pion contribution is 

{P2)p{pi)a{pt +pf^) = 2 [Pp + I) (P„ - y) = 2PpP^P^ + 0{q), 

and the "best" structure (not changing when q varies) is PaPpP^ ■ The simplest 
way to extract the relevant amplitude (hereafter referred to as T) is to multiply 
T^p{pi^P2) by nfj^n^n^, where n is a lightlike vector orthogonal to q. The 
property = kills all the structures containing gai3 and other y-terms, 
while the property (nq) — kills all the structures containing any g-factor. 
Thus, the amplitude we are going to analyze is 

nplplq')-^I^T^p{pi,P2)- (191) 

7.2 Dispersion representation 

To extract information about the form factors of physical states, one should 
incorporate the double dispersion relation 

2 2 2n 1 f°° ^ p{si,S2,q^) „ 

1 (p,,Pn,q ) — ^ / d,si / as2 7 ^77 ^+ subtractions . 

Jo Jo [si - Pi){s2 - Pi) 

(192) 

The subtraction terms, in general, are polynomial in pi and/or p^- They dis- 
appear after one applies B{p^ M^)B{p\ — > Ajf|), the Borel transformation 
in both variables 0. The double Borel transform $ = B1B2T is related to the 
double spectral density p{si,S2,q^ = —Q^) by 

(193) 

7.3 Operator product expansion 

In contrast to the two-point correlators Il{Q^) analyzed in the previous Section, 
the three-point amplitudes T(p^,p2i'?^) have three arguments. To incorporate 

In our specific case, the subtractions are exhausted by a constant, and a single differenti- 
ation with respect to any variable is sufficient to remove it. 
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a) b) c) d) 



Figure 8: Structure of OPE in case of symmetric kinematics. 

asymptotic freedom, one should take the invariants p\,p\ corresponding to the 
"hadronized" channels to be negative and sufBciently large, say, IpiMpII > 
1 GeV^ . Depending on the value of the third invariant that works as an 
external parameter (it also should be negative, of course), there are three 
essentially different situations. 

• The simplest is the symmetric case when is of the same order of magnitude 
as \pI\ApI\- 

\pl\-\pl\-Q'-^^'■ 

The operator product expansion for T in this case has the following structure: 

T{pIpI q') = E AwC. f 4. 4> 4) m (194) 

where Ci/{iJ?Y^ is a perturbatively calculable short-distance coefficient func- 
tion, {Oi) is the vacuum matrix element of a local operator , i.e., some con- 
densate term and di is the dimension of Oi in mass units (see Figj|). 

• The case when is small is more complicated: 

\pl\-\pl\-^^^ Q'«t^'. 

The coefficient functions Ci calculated in the symmetric kinematics, might 
become singular in this limit because of {fi^ /Q^)^- or log(/.t^/(3^)-terms. The 
operator product expansion in such a situation has a modified form: 
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where the coefficient functions C, are regular in the — > Umit, Ilk{Q^) is 
the correlator of the electromagnetic current J and a local operator Ok 

dk is the dimension of Ok in mass units. The coefficient functions Ek arc (the 
Fourier transforms of) those for the operator product expansion of two axial 
currents 

T{jj+)=Y,EkOk. 

k 

In the modified OPE for the three-point function, the contributions of the 
first type correspond to a situation when large momentum flow connects all 
three external vertices, while those of the second type describe the situation 
when large momentum flows only between the vertices x and y related to the 
ax;ial currents. The momenta in the J-channel are of the order of Q, i.e., 
small, and the correlator Ilj{Q'^) is a nonperturbative object accumulating 
long-distance information. 
• In the opposite limit of large 

a simple expansion in powers of |, l/l^i (*-e-)in might be destabi- 

lized by large ((5^//u^)"-terms. Appearance of terms growing with is a bit 
surprizing, since form factors should decrease as oc'. But the appearance 

of such terms is just the artifact of the 1/ fi^ expansion: even a very reasonable 
function like 1/{Q^ + fji^), vanishing as oo, produces a divergent series 

when expanded in Thus, a possible way out in the large-Q^ limit is a 

resummation of the parametrically enhanced contributions of (<5^/At^)""t;ype. 
In fact, it can be established, that such contributions result from the Taylor 
expansion of the original nonlocal combinations (!^(0)f^(a;)) (nonlocal conden- 
sates), and the idea is to construct sum rules without expanding the nonlocal 
condensates, treating them as some phenomenological functions describing the 
propagation of quark and gluons in the physical vacuum. 

7.4 Perturbative contribution 

The structure of the power corrections in the operator product expansion, as we 
discussed above, depends on the particular Q^-region under study. However, 
the starting point of any expansion is the perturbative contribution, in our 
case corresponding to the triangle diagram plus radiative corrections to it. 
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Using Feypman parametrization, one can rather easily obtain the following 
representatioiO for the double Borel transform of the lowest-order perturbative 
triangle diagram (see Fig.||a) 

"''P {~ {1 - xYMi + hq) ] ^^1' '^'^ 

where I{x, , ,111^, Cq) is the factor containing the dependence on the 
quark masses 

I{x, Ml Ml ml, e,) = e„ exp |- + (^_L + | + ^ d}. 

(197) 

In view of extreme smallness of the u- and d-quark masses, this factor can be 
safely set to unity. The x-variable may be interpreted as the fraction of the 
total momentum P carried by the passive quark (active is the quark interacting 
with the EM current). This integral representation is very convenient to study 
the behavior of $P<=''*(Mf , M|, Q^) for smaU and large Q^. in the low-Q^ limit 
we get 

*-'(M?, Ml, 0=) ^ ^,^, {1 - (198) 

+ ^-*(iV + l))}. 

The first two terms of the expansion can be obtained by simply expanding the 
exponential factor in under the integral sign. However, the third term of 
such an expansion produces a logarithmically divergent integral 

f' If Q' y , 

■■■Jo 2! [Mi + MiJ l^x"'"' 

and higher terms have power divergences when integrated in the region x ^ 1. 
This region corresponds to a situation when the whole large momentum P 
is carried by the passive quark, and the active quark carries small momen- 
tum proportional to q. This is precisely the situation producing additional 
terms in the operator product expansion at small Q^. The logarithmic factor 
log((5^), singular when 0, just signalizes that there are long-distance 



E 

AT n 



Af2 



iV! 



log 



Ml 



Mi 
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effects (quarks propagate over distances ~ l/Q) limiting the applicability of 

perturbation theory. 

In the large-Q^ limit one should expand ^P''''^{Mf, M|, Q^) in l/Q^ rather 
than in Q^: 



+ \ E (-)"(iV + l)(iV + 3)! (^M±M)''|. (199) 

Thus, the perturbative contribution vanishes like l/Q"* as — » oo. Combining 

the spectral representation for <I>(M^, M|, Q^) and the explicit expression for 
it, we get the equation for the perturbative spectral density p^"^''*(si, S2, Q^) 

1 [°° dsi f°° ds2 , Si S2 t 2n 

Solving it (this amounts to a double inverse Laplace transformation) gives the 
perturbative spectral density: 



2\3 



2! 3! 



V(si+S2 + (32)2_4siS2 



(201) 



In the low-Q^ limit the perturbative spectral density can be expanded in Q'^: 

p^^'^*(si,s2,q')|q2^o = - S2)e{s,)e{s2) + 

+ ^(si + S2)<5"(si - S2)9{sr)e{s2) + .... (202) 
For large one can expand it in 1/Q^: 

(Sl, S2, g )|q2^<^ = - 6l(Si) e{S2) < 4 + . . 

(203) 

Thus, for small Q^, the perturbative spectral density is concentrated on the 
line Sl = S2 ■ there are no transitions between states with different masses. 
Along this line, the spectral density does not vary: it has no dependence on 
Sl + S2- In the opposite limit, when oo, the situation is, in a sense, 
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reversed: the density is nonzero in the whole region si > 0,S2 > 0, there 

is no dependence on (si — S2), and the density grows like si + S2. The last 
observation means that higher states are relatively more important for large 
Q"^ than for the small ones. 

7.5 Sum rule for the pion form factor at moderate 

The physical spectral density p{si,S2,Q^), of course, differs from its pertur- 
bative analog. It has a rather complicated structure on the si,S2-plane. In 
particular, p(si, S2, Q^) contains the term corresponding to the pion form fac- 
tor 

P^^{si,S2, Q^) = Tr^f^F^{Q^)6{si - ml)6{s2 - ml). (204) 

In addition, it contains the contributions corresponding to transitions between 

the pion and higher resonances, and also the terms related to elastic and tran- 
sition form factors of the higher resonances. Constructing the two-dimensional 
(si,S2) analog of the simplest ansatz "lowest state plus continuum" we can 

treat all the contributions, except for the Ptttt, as "contimmm", i.e., assume 
that p{si,S2,Q^) = p^'^''*(si, S2, Q^) outside the square (0, so) x (0,,so): 

P(S1, S2, Q^) = Pnnisi,S2, Q^) + (1 - 6{S1 < So)e{S2 < So)) p^"=''*(si, S2, Q^). 

(205) 

Now, transferring the continuum contribution to the right-hand side of the sum 
rule and using the operator product expansion in the symmetric kinematics 
{Q^ ^ Ml, M|), we obtain the QCD sum rule for the pion form factor: 



/X(<3') = ^^"'dsi^'°ds2exp(-^i^ ) pP-*(si,S2,Q') + 



To treat initial and final states on equal footing and to simplify the ana 
we took = Ml . Unfortunately, the double integral in this sum rule cannot 
be calculated explicitly for arbitrary Q"^. However, it is instructive to consider 
the formal limit = 0: 

fFm-^^il- ^ + "-^^^^ ^ (207) 

J^i'-W - 4^2 2 r J + 127rM2 +81 M4 • ^^^'> 

This should be compared to the sum rule: 

f2 _ ™' ^-so/m^\ : (^s{GG) 176 7ra,(gg)2 

^ Att^V ^ J + 127rm2 + 81 ' ^^"^^ 



53 



The condensate terms look almost identical, but comparing the perturbative 
terms we observe that M^, the Borel parameter of the form factor sum rule, 
should be larger by factor 2 than m^, the Borel parameter of the sum rule: 
— 2m?. After this change, the (GG')-term in the form factor sum rule is 
multiplied by ^ and the (^g)^-term is multiplied by -j 

Thus, we obtained a sum rule, which looks like that for /^, but with the gluonic 
condensate term smaller by factor 2 « (1.4)^, and the quark condensate term 
smaller by factor ~ 3.5 « (l-5)'^. Since it is the values of the condensates that 
(after the fitting procedure) determine all the hadronic scales, we conclude that 
all the hadronic paramerters having the dimension of {mass)'^, i.e., sq and the 
combination f^F.^{Q), extracted from this sum rule are by a factor 1.4 — 1.5 
smaller than sj and respectively]^. To get sq closer to sj, one should restore 
the matching between the perturbative and the condensate terms, i.e., decrease 
the perturbative term by a factor 2 to 4. This allows us to make an educated 
guess that QCD sum rule for the pion form factor will most effectively work 
in the intermediate region where the form factor varies between 0.5 and 0.25, 
i.e., for Q2 between 0.5 GeV^ and 2GeV^ 

The readers can check these guesses by explicit fitting procedure. To this 
end, one should plot the combination f'^FT^{Q'^) as a function of the Borel pa- 
rameter for different values of the effective threshold s^. As a "true" value 
of So one can take the value for which the curve is constant as Af ^ oo. One 
would observe that the value of sq obtained in this way slowly grows with Q^, 
from 0.6 GeV^ for = 0.5 GeV^ to 1.0 GeV^ for = 3 GeV^ The growth of 
So reflects the fact that, when the value increases, the condensate contribu- 
tions remain constant (with the quark term even slowly growing) whereas the 
perturbative term decreases. The constancy of the condensate contributions is 
an artifact of our approximation. If one includes operators of higher dimensions 
in the OPE, there appear terms diminishing the total condensate contribution 
for large Q^. These terms, however, have the structure (Q^/M^)", and one 
should resum them to get a reasonable result. However, if one restricts the 
analysis to the region 0.5 GeV^ < Q2 < 3GeV^ the predictions of our sum 
rule agree with existing pion form factor data within 10% — 20% accuracy. 

"^One should not worry about the fact that a straightforward extrapolation of the form factor 
sum rule to the point Q-^ = 0, gives F,r(0) ~ 0.7. This sum rule is valid only for sufficiently 
high Q^, and, as we emphasized, there appear additional terms in the OPE which will 
increase the condensate contributions to make them equal to those in the sum rule and 
recover the -F7r(0) = 1 normalization condition. 
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7. 6 Local quark-hadron duality for the pion form factor 

The approach described above allows to obtain the form factor at different 
values of by independent fitting procedure, point by point. But it is 
tempting to see the Q^-curve as a whole. This can be done by looking at 
the F^(Q^)-curves at fixed sq for different values of the Borel parameter M^. 
For M2 > 1.5GeV^ all the curves are close to each other, approaching the lim- 
iting (sQ-dependent) form as oo. The nature of this convergence is very 
simple: the condensate terms disappear from the sum rule in the Af^ — + oo 
limit, and the limiting curve corresponds to the local duality relation: 

flF^{Q-^) = — dsi pP"'-*(si,S2,Q')ds2. (210) 

Jo Jo 

Note, that the function pP'^^*{si, S2,Q^) describes the transition from a free- 
quark (qq) state with invariant mass si to a similar state with mass S2- The 
local duality relation states, essentially, that one can calculate the pion form 
factor by averaging the form factors of such transitions over the appropriate 
duality interval sq. In the two-point case, the local duality produces the rela- 
tion between the pion decay constant and the duality interval sq 

So = 47r2/J ~ 0.7GeV^ (211) 

We recall that the duality interval is the effective threshold for production of 
higher resonances, and is very natural that it is just in the middle between 
ml wJland « 1.6 GeV^. Using the explicit form for /9P'^''*(si, S2, Q^), we 
obtaiiEil 

p ( r)2\\local duality _ i 1 -|- 6so/ Q (0^0\ 

This formula gives a correct value for the form factor at = 0, but the 
predicted slope is wrong: 

p(Q^\\localduaUty ^ ^ _ 3 _Q_ ^ ^ ^ ^ .2I3) 

4 V«o 

This is just the same situation we encountered applying the local duality to 
the form factor of the lowest oscillator state. The reason is that there is a 
serious mismatch between perturbative and physical spectral densities: the 
perturbative density is concentrated in the narrow stripe |si — S2I < Q near 
the si = S2 line, while physical density, in addition to the elastic form factors 
(points at the si = S2 line), contains also transition form factors (line si = 
0, S2 > s(3'^*'"''2«^o''')^ etc.) corresponding to the regions very far from the si = 
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■32 line. In other words, our model for the hadronic spectrum is too rough at 
small nonzero Q^. The perturbative density in this region can be approximated 

by 

/^'■*(si,S2,Q2)^^0(|5i-S2|<Q). (214) 

Integrating it over any line si + S2 =const., one obtains a result analytic in 
Q^. In our model, however, there are small triangles not included into the 
integration region. Their area is ~ Q^, and the density is ~ l/Q, so the 
resulting contribution (to be subtracted from 1) is 0{Q) = 0{y/Q^). To avoid 
this, one can try a "triangle" model for the higher states: "pion + (continuum 
outside the triangle forrnnd by the hne (si + S2 = Sq)" . The local duality 
relation in this case giveO 

Tp I /W2\trianglel.d. /'91 

^^^^ " 8^2;2(l + Q2/25„)2- ^^^^i 

To get the correct normalization -^^(0) = 1 one should take Sq — Sn'^f^ = 2so. 
However, the slope in this case is too small -F^'' ''^ (Q^)|so=2so = l^Q^/2so, by 
factor 2 smaller than the experimental one, and the curve goes well above both 
the experimental points and the curve corresponding to the "squared" form of 
the local duality. To match with the latter, one should take 5*0 = V^sq. In 
this case the area of the local duality region is the same. Outside the small-Q^ 
region, there is good agreement between the two local duality curves and the 
data approximated by the fit Ff'\Q^) «!/(! + Q^/0A7). 



7.7 One-gluon exchange contribution 

The contribution into the pion form factor we calculated using the lowest order 
0{oP^) term in the perturbative spectral density corresponds to the so-called 
"soft diagram". Since this contribution decreases as l/Q^ for large Q^, it is 
normally ignored in the perturbative QCD analysis concentrated on the study 
of the "hard gluon exchange diagram" that has the asymptotic behavior 

dictated by the quark counting rules. 

However, the soft contribution as given, e.g., by the local quark-hadron 
duality, is very close to the experimental data leaving not much place for any 
other contribution. To estimate the contribution due to the "hard gluon" ex- 
change diagrams, one should include the 0(Q!s((5^)-terms in the perturbative 
spectral density. To simplify the analysis, one can use the local duality approx- 
imation. Still, the two-loop calculation is rather complicated, but the results, 
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a) 



b) 



Figure 9: Gauge-invariant set of {qq)^ diagrams. 



with a rather high accuracy, can be approximated by a simple interpolation 



between the = value (related by the Ward identity to the 0{as) term of 
the 2-point correlator) and the asymptotic behavior F^^{Q'^) = 8Tras{Q^)f^/Q'^. 

7.8 Pion form factor at small 

As we discussed earlier, if the QCD sum rule derived in the moderate-Q^ 
region is simple-mindedly extrapolated to the point = 0, one obtains the 
expression 



The condensate terms have the eoefficients essentially smaller than it is nec- 
essary to enjoy the simple Ward identity prediction F.^{Q'^) — 1. The latter 
requires that the QCD sum rule for /^iv(0) should coincide with that for 
The Ward identity, being just the statement that the total electric charge of 
the pion is equal to that of the electron, should not be violated. Hence there 
should exist the additional condensate contributions "visible" at = (to 
satisfy the Ward identity) and dying rapidly for > (to reproduce the 
intermediate-Q^ sum rule). 

To get a feeling about how these terms might look like, consider the ag (qq)'^ 
diagrams. To enjoy the consequences of the Ward identity, one should have a 




(216) 




(217) 



that differs from the sum rule for 




(218) 
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gauge invariant set of diagrams. This can be achieved by inserting the photon 
vertex into all possible lines, including also the "external" lines corresponding 
to quarks going into vacuum. However, in perturbation theory, this brings 
in terms ^ singular in the — > limit. Such a "coefficient function" 

cannot be calculated perturbatively. Thus, one should separate the large- and 
small-momentum parts of the diagram 

T{p',q^)^C{p')®n{q^) (219) 

and treat the Q^-dependent part as a correlator of the original electromagnetic 
current J and some local operator constructed from the quark fields going "into 
vacuum" (see Fig.^): 

n(g2)^ j e''i^T{J^{Q)0{x)))d\. (220) 

The question now is whether it is possible to calculate contributions like the 
correlators 11 (g^) at small ql There are essentially two types of these correlators 
depending on the type of the operator O. First, one should take into account 
that under the T-product sign, one cannot throw out the operators apparently 
vanishing due to the equations of motion, e.g., those containing^ acting on 
the quark field, since there appear contact terms. In this case n(g^) is a 
constant proportional to the quark condensate {qq) . In all other cases one can 
(approximately) calculate n(g^) using the following QCD sum rule strategy: 
• Calculate n(q^) for large using the operator product expansion. This 
normally gives 

n(,^ = -Q^)^M + fc<fg^ + .... 



• Use the dispersion relation 

p{s) ds 



1 

m') - - / 

I" Jo 



• Assume an appropriate ansatz for the spectral density: 

p{s) — nAS{s — rn?p) + nBS{s — rnfj). 

Since the perturbative density is zero in this case, the higher states (lying 
above the p-meson) are approximated by an effective resonance R. 

• Requiring the best agreement between the two representations for n(Q^), 
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extract the values of the constants A, B and 



The result is 



n(Q2 



A 



B 



and one can use it for small Q^. 

In this way one can obtain the QCD sum rule for the pion form factor 
valid in the low-Q^ regioiH: 



/^-F!n-((3^) — ^^perturbativeparf^ 



16 



nas{qq) 



81 



4 + 6 



1.6 m2 



asjGG) 



0.6 to| 



2{ml 



(221) 



This sum rule possesses all the required properties: 

• For — Q its left hand side exactly coincides with that of the sum rule. 

• additional terms die out when gets large. 

• The curve for the pion form factor calculated with the help of this sum rule 
in the region of small , at the matching point — rn^, agrees within 10% 
with the curve obtained from the moderate-Q^ sum rule. 

• the pion charge radius, extracted from the low-Q^ sum rule, 



,2n1/2 



5_R = 0.66 ± 0.03 fm 



is in good agreement with the experimental value 



.2 \ 1/2 I 



-,p. = 0.636 ± 0.036 /m. 



8 Conclusions 

Thus, the QCD sum rules give an accurate description of the pion form fac- 
tor in the whole region where the data exist. They allow one to trace its 
behavior from the normalization point = 0, through the low-Q^ region, 
to the region of moderately large Q^. Everywhere the QCD sum rule results 
are in agreement with experiment within the accuracy of the method. The 
QCD sum rules unambigously demonstrate that, in the experimentally acces- 
sible region, the form factor is dominated by the soft contribution, which is 
not calculable within a purely perturbative approach. They show that the 
one-gluon-exchange mechanism, though dominant in the asymptotic limit, is 
of minor importance in the region Q'^ < 4 GeV^ and, maybe, till even larger 
values of Q^. 
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